Self-trapping of Bose-Einstein condensates in optical lattices 
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The self-trapping phenomenon of Bose-Einstein condensates (BECs) in optical lattices is studied 
extensively by numerically solving the Gross-Pitaevskii equation. Our numerical results not only 
reproduce the phenomenon that was observed in a recent experiment [Anker et al, Phys. Rev. Lett. 
94 (2005)020403], but also find that the self-trapping breaks down at long evolution times, that is, 
the self-trapping in optical lattices is only temporary. The analysis of our numerical results shows 
that the self-trapping in optical lattices is related to the self-trapping of BECs in a double-well 
potential. A possible mechanism of the formation of steep edges in the wave packet evolution is 
explored in terms of the dynamics of relative phases between neighboring wells. 

PACS numbers: 03.75.Lm,03.75.Kk,05.45.-a 
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I. INTRODUCTION 



Progress in recent years has shown that a Bose- 
Einstein condensate(BEC) in an optical lattice is a fas- 
cinating periodic system, where the physics can be as 
rich as in fermionic periodic systems, the main sub- 
ject of condensed-matter physics. In such a bosonic 
system, people have observed well-known and long pre- 
dicted phenomena, such as Bloch oscillations [lj and the 
quantum phase transition between superfluid and Mott- 
insulator[2j. More importantly, there are new phenom- 
ena that have been either observed or predicted in this 
system, for example, nonlinear Landau-Zener tunneling 
between Bloch bands 0, and the strongly inhibited 
transport of one dimensional BEC in an optical lattice^. 

Another intriguing phenomenon, self-trapping, was re- 
cently observed experimentally in this system|6(. In this 
experiment, a BEC with repulsive interaction was first 
prepared in a dipole trap. By adiabatically ramping up 
an optical lattice, the BEC was essentially transformed 
into a Bloch state at the center of the Brillouin zone. 
With the optical lattice always on, the BEC was then 
released into a trap that serves as a one dimensional 
waveguide. The evolution of the BEC cloud inside the 
combined potential was studied by taking absorption im- 
ages. When the number of atoms in the BEC is small, 
say around 2000, the BEC wave packet was found to ex- 
pand continuously (which is expected). However, when 
the number of atoms was increased to about 5000, it 
was observed that the BEC cloud stops to expand af- 
ter initially expanding for about 35ms (see Fig^). This 
is quite counter-intuitive. Without interaction, a wave 
packet with a narrow distribution in the Brillouin zone 
expands continuously inside a periodic potential. One 
would certainly expect that with a repulsive interaction 
between atoms the BEC cloud spread faster. This exper- 
iment showed the contrary: if the cloud is dense enough, 
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Figure 1: The width of the BEC wave packet as a function of 
time for N = 2000 and N = 5000. N is the number of atoms 
in the BEC. The solid lines are our numerical results while 
the circles and squares are experimental data from Ref.Q. 



it self-traps and stops spreading. 

To understand this intriguing phenomenon, we have 
carried out extensive numerical study of this system with 
the one-dimensional(lD) Gross-Pitaevskii equation. Our 
results match quite well with the experimental data as 
shown in Fig^ When the atom number N in the BEC 
is 2000, the agreement between our numerical results and 
the experiment is excellent; when N = 5000, our re- 
sults are about 40% larger than the experimental data. 
The discrepancy in the latter case is likely caused by the 
higher density: with higher density the lateral motion of 
the BEC cloud may become more relevant to the longi- 
tudinal expansion; however, the lateral motion is com- 
pletely ignored in our numerical study as we use the ID 
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Gross-Pitaevskii equation. 
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Figure 2: The width of the BEC wave packet as a function of 
time for N = 1000, 2000, 3000, 4000, 5000, 7000, and 50000. 

Very interestingly, we find that the self-trapping is 
temporary. After a sufficiently long evolution time, the 
self-trapping breaks down and the wave packet starts to 
expand again as seen in Fig|2 Since the break-down time 
is much longer than the observation used in the current 
experiment @, these results need to be verified in future 
experiments. This breakdown of self-trapping is likely 
caused by the leakage of atoms at the outmost wells, or 
as we shall call it the dripping effect. Furthermore, our 
numerical results show that the steep edges do not nec- 
essarily lead to the self-trapping and always appear in 
the wave packet evolution, independent of its denseness 
of the BEC cloud. This is different from Kef. (i . where 
it was pointed out that the steep edges appearing at the 
two sides of the wave packet are crucial for the appear- 
ance of the self-trapping. The wave packet evolution in 
the quasi-momentum space is also studied. We find that 
the wave packet localizes largely near the center of the 
Brillouin zone; no major interesting features can be iden- 
tified during the evolution. 

Besides, we have analyzed our numerical results in de- 
tail, in particular, in terms of the relative atom num- 
ber difference and the relative phase difference between 
neighboring wells. With such an analysis, we have con- 
firmed the previous study that the self-trapping ob- 
served in optical lattices is closely related to the self- 
trapping of a BEC in the double- well potential 0, El • 
From this analysis, we have also explored a possible 
mechanism of the formation of steep edges. 

Our paper is organized as follows. We shall first 
present a brief description of our numerical method. We 
then describe how the wave packets evolve in our numer- 



ical simulation. Afterwards, we analyze our numerical 
results in an attempt to understand our numerical re- 
sults. The analysis is done from the angle of a BEC in a 
double- well potential. Summary and some discussion are 
given at the end. 



II. NUMERICAL METHOD 

To model the experiment, we use the following Gross- 
Pitaevskii(GP) equation 

d h 2 
ih—il>(r, t) = -—V 2 ip(r, t) + V Q cos{2k L x)ip(r, t) + 
at 2m 

+V wg (r)xP(r,t) + —^\^(r,t)\ 2 ^(r,t) , (1) 

where m is the atomic mass, a s is the s-wave scattering 
length, is the wavelength of the laser that generates 
the optical lattice, and V wg (r) describes the waveguide 
potential. Due to the tight confinement perpendicular 
to the optical lattice from the waveguide potential, the 
dynamics of this system is largely one-dimensional. This 
allows us to integrate out the two perpendicular direc- 
tions and reduce the above GP equation to 



.dip(x, t) 
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V 2 ip{x, t) + V cos(x)ip(x, t) + 
ux 2 i>(x,t) + g\ip(x,t)\ 2 ip(x,t) , (2) 



were we have made the equation dimensionless. In doing 
so, we have x in units of 1/2/cl and t in ra/Ahk\. The 
strength of the optical lattice is given V — Vo/16E r with 
E r = h 2 k\/2m being the recoil energy. For the nonlinear 
interaction, we have 
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where N the total number of the BEC in the harmonic 
trap and lu± is the transverse trapping frequency of the 
waveguide. The other frequency u> is so chosen that 
the initial rms-width of BEC wave packet is 7.6/Ltm as 
in the experiments 0- This width corresponds to about 
100 wells occupied. The wave function ip is normalized 
to one. In our numerical simulation, we use the follow- 
ing values from the experiment 0, A = 2ir jk^ — 783nm, 
V = WE r , and uj±_ = 2n x 230Hz. 

To simulate the experiment, we prepare our initial 
wave function ip to be the ground state in the combined 
potential of V cos(x) + \u>x 2 . This is achieved by inte- 
grating Eq.J2J with imaginary time. In the experiment, 
the waveguide potential also has a longitudinal trapping 
frequency at u>\\ = 27rHz, which is very weak and can be 
ignored. Therefore, after obtaining the initial wave func- 
tion, we completely remove the longitudinal trapping and 
let the wave function evolve according to the following 
equation 

.dip 1 



at 



V 2 ip + V cos(x)tp + g\-ip\ z ip ■ 



(4) 
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The evolutions are subsequently recorded and analyzed. 



III. WAVE PACKET EVOLUTION 

With the above method, we have computed the evolu- 
tion of the wave packets for different numbers of atoms in 
the BEC. As indicated in Eq.©, the number of atoms in 
the BEC translates into the nonlinear parameter: larger 
the atom number N stronger the nonlincarity (or the 
repulsive interaction). Fig|2 illustrates how the width 
of a wave packet evolves for different atom numbers. 
It is clear from this figure that, when the BEC is di- 
lute and has small atom numbers, N < 2000, the wave 
packet expands continuously without stopping as one 
may have expected. Also as expected, in this range, 
when the number of atoms increases, the expansion be- 
comes faster. The evolution becomes very different when 
the BEC is denser. We see in FigElthat for N = 3000, 
the wave packet expansion slows down around 70ms and 
becomes slower than the wave packet for TV = 2000 
around 85ms. This means that we have slower expan- 
sion for a denser cloud of repulsive interaction, a rather 
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Figure 3: Time evolution of the wave packet density for N 
5000. 



counter- intuitive result. As the cloud gets denser with 
more atoms, the expansion slows down further. Around 
N = 5000, there even appears a plateau where the cloud 
stops expanding and becomes self-trapped as observed in 
the experiment |(| . Fig[2 illustrates a key point in this 
intriguing self-trapping phenomenon: it does not happen 
in a sudden and it is a gradual process. Before it hap- 
pens, the wave packet expansion already slows down for 
higher enough densities. 

What is more interesting is that, in our numerical 
simulation, the wave packet continues to expand after 
pausing for 30-40ms. For N = 5000, the expansion 
re-starts at ~ 85ms, just beyond the longest observa- 
tion time in Ref. 6]. Therefore, this continued expansion 
awaits for verification in future experiments. Neverthe- 
less, the counter-intuitive phenomenon, denser clouds ex- 
pand slower, persists even after the expansion re-starts 
as we can see in Fig|2 In the next section, we shall of- 
fer an explanation of this self-trapping phenomenon and 
explain why it is only temporary. Just out of curiosity, 
we also computed the case of very large atom number 
N = 50, 000; we find that the wave packet is almost 
never seen to spread. This is likely due to that the self- 
trapping lasts too long to be observed in our numerical 
simulations. 

Shown in Fig|3|are some snapshots of the time evolu- 
tion of the wave packet for N = 5000. Around t = 30ms, 
steep edges are seen growing pronounced on both sides 
of the wave packet. Moreover, in the subsequent evolu- 
tion, the positions of the steep edges do not move out any 
further. When compared with Fig^ it is clear that this 
appearance of the steep edges coincides with the non- 
spreading of the wave packet. Therefore, it seems that 
the steep edges seen in Fig |3] signal the emergence of the 
self-trapping as suggested in Ref. 6]. The following re- 
sults indicate otherwise. 

Fig^ shows the time evolution of the wave packet for 
N = 2000. Before 80ms (about the longest experimental 
observation time in Ref. 6]), there are no steep edges. 
However, around 85ms, the steep edges begin to ap- 
pear and grow more and more pronounced as the evo- 
lution goes on. What is different from N — 5000 is that 
these steep edges continue to move out during the time 
evolution and the width of the wave packet also grows 
with time. There is no self-trapping. This clearly shows 
that the steep edges do not necessarily lead to the self- 
trapping of the wave packet. 

We have also computed how the wave packet evolves in 
the quasi-momentum space. To achieve this, we expand 
the wave packet in terms of the Bloch waves belonging 
to the lowest Bloch band of the linear system with the 
periodic potential cos a;. The results are plotted in FigHj] 
where we do not see a large population around k = 1/4. 
Therefore, the link between the formation of steep edges 
and the population at k = 1/4 believed in Ref. [a] is not 
established here. 

Another quantity that can be used to characterize the 
self-trapping phenomenon is the nonlinear energy of the 
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BEC, which is given by J ^\ip(x, t)\ A dx. If the wave 
packet expands continuously, the nonlinear energy should 
decrease with the expansion. If there is self-trapping, i.e., 
the wave packet stops to grow, then the nonlinear energy 
should remain largely constant. Indeed, this is the case as 
shown in Fig0 where the nonlinear energy docs not de- 
crease and only fluctuates slightly when the self-trapping 
occurs. 



IV. RELATIONSHIP TO THE SELF-TRAPPING 
IN A DOUBLE- WELL 



It has been known for a while that the self-trapping 
also occurs for a BEC in a double- well potential 0, H, Hi 
ITfij . It is then natural to ask whether the self-trappings 
in these two different systems are related to each other. 
Our analysis shows that these two are closely related. 
This was already noticed in Ref.Q based on numeri- 
cal results with a tight-binding approximation of the GP 
equation; here we offer a more detailed analysis with nu- 
merical results with the full GP equation. Our analysis 
leads to a possible explanation why the self-trapping ob- 
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Figure 5: Time evolution of the wave packet density for N = 
2000 after t = 80ms. 



served in the optical lattice is only temporary and how 
the steep edges form. For the sake of self-containment 
and introducing new parameters, we briefly review the 
self-trapping in the double-well situation. 

The Hamiltonian governing the dynamics of a BEC in 
a double-well potential can be written as[~ 



H c la 



ssical 
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v\/l — s 2 cos 6, 



(5) 



Figure 4: Time evolution of the wave packet density for 
2000 before t = 80ms. 



where 9 = 9b — a is the relative phase between the two 
wells a and b while s is the fractional population dif- 
ference s = (N b - N a )/(N a + N b ) with N a and N b be- 
ing the number of atoms in wells a and 6, respectively. 
Previous studies 0, IE show that there are two types 
of self-trapping in this system, depending on the ratio 
£ = c/v and the population difference s. They are: (1) If 
1 < ^ < 2 and the relative phase 9 is around 7r, then the 
self-trapping occurs when s > 0.5. This is called 'oscilla- 
tion type' self-trapping. (2) If £ > 2 and s > 0.5, another 
type of self-trapping emerges with the relative phase 9 be- 
tween the two wells increasing with time. Therefore, it 
is called 'running phase type' self-trapping. 

Any pair of neighboring wells in the optical lattice can 
be viewed as a double-well. To establish the link between 
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Figure 6: Time evolution of the wave packet in the quasi- 
momcntum space for TV = 5000. 



of these four pairs, marked by A and B in Fig [SI are 
located right at the two edges of the wave packet. (The 
other two are just nearby, for clarity we do not mark 
them.) It seems to suggest that these two self-trapped 
pairs of double-wells serve as two dams stopping the flow 
of atoms to the outside. Therefore, the self-trapping in 
the optical lattice appears just alternative manifestation 
of the self-trapping in the double-well system. 

To firmly establish such a link, we have also examined 
the case TV = 2000, where there is no self-trapping. Fig|5| 
shows the values s and £ for the wave packet with TV = 
2000 at t = 110ms. This is the time when the steep edges 
have already developed. We see from the figure that all 
the values of £ are smaller than one: the self-trapping 
conditions of the double-well system are not satisfied by 
any pair of neighboring wells in the optical lattice. 

The above analysis leads to the following conclusion: 
when there is self-trapping in the optical lattice, there are 
neighboring wells that satisfy the self-trapping condition 
of the double-well system; when there is no self-trapping 
in the optical lattice, any pair of the neighboring wells in 
the lattice does not satisfy the double-well self-trapping 
condition. So established is a solid link between these 
two self-trapping phenomena. One intuitive way of un- 
derstanding of this link is such. Once the self-trapping 



the self-trappings in the optical lattice and the double- 
well, we need to compute £ and s for each pair of the 
neighboring wells in the optical lattice for a given wave 
packet. The details of how £ is computed for neighboring 
wells in an optical lattice can be found in Appendix A. 

Fig [5] shows one set of such calculations for a wave 
packet with TV = 5000 at t = 40ms, which is the time 
when the self-trapping happens. It is clear from FigQJa) 
that there are four pairs of double-wells whose £ and s 
satisfy the condition for the "running phase type" self- 
trapping in the double-well system. Furthermore, two 
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Figure 7: The dynamic of nonlinear energy of each atom for 
TV = 2000 (dashed) and TV = 5000 (solid). 
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Figure 8: The values of £ and s of a wave packet at t = 40ms 
for TV = 5000. (a) The values of £ at different wells (The £ at 
the nth well is for the double wells composed of the nth and 
n + 1th wells, (b) Phase diagram (or distribution) of £ and s 
for this wave packet. The self-trapping occurs at t — 40ms. 
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happens in some pairs of neighboring wells around the 
edges of a wave packet, these self-trapped double-wells, 
behaving like "dams" , stop the tunneling of atoms to- 
wards outside, causing the non-spreading of the wave 
packet. In Fig^]we have plotted a series of "phase dia- 
grams" , where the distribution of the s-£ pairs is shown. 
The square region bounded by the dashed line in each 
panel is the area where the self-trapping conditions are 
satisfied. We can see clearly from this figure that the 
self-trapping happens from about t = 40ms to 80ms for 
N = 5000. We have also checked the cases of N = 7000 
and N = 50, 000 and reached the same conclusion. 

This link not only explains why the self-trapping oc- 
curs in the optical lattice but also offers a possible mech- 
anism why the self-trapping is temporarily lived. At the 
outmost wells, the density of the BEC is very low and 
the self-trapping conditions of the double- well system can 
never be satisfied. As a result, the atoms will tunnel to- 
wards outside. The amount of atoms tunneling out is 
very small and has not much effect on the evolution of 
the whole cloud. However, for long evolution times, this 
small amount of "dripping" can lead to the significant 
decreasing of atom numbers in the wells, thus destroy 
the self-trapping. This is similar to that small cracks can 
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Figure 10: Phase diagrams (or distributions) of s and £ for 
wave packets at different evolution times. N = 5000. 
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cause the collapse of a dam in a long time. 

Based on the link between these two self-trappings, it is 
also possible to understand why the BEC cloud expands 
slower around N = 3000 than N = 2000 seen in Fig|2 
As one can imagine, when the cloud density increases, 
some pairs of the wells will get close to satisfy these self- 
trapping conditions and eventually satisfy them. For the 
medium densities, e.g., N — 3000 there should be a few 
pairs of the wells that satisfy the conditions just barely. 
As a result, the self-trapping conditions can be easily 
or quickly destroyed by the "dripping" effect mentioned 
above. However, as the cloud expands, the self-trapping 
conditions can again be satisfied by some pairs of well 
further inside and then destroyed again. This on-and- 
off process can dramatically lead to slowing down of the 
cloud expansion. What is a pity is that this straightfor- 
ward picture is hard to be corrobarated by our numerical 
computation because the values s and £ for neighboring 
wells can only be computed approximately. Alternative 
methods may be needed to verify this picture. 



Figure 9: The values of £ and s of a wave packet at t = 110ms 
for N = 2000. (a) The values of £ at different wells (The £ at 
the nth well is for the double wells composed of the nth and 
n + 1th wells, (b) Phase diagram (or distribution) of £ and s 
for this wave packet. The time 110ms is when the steep edges 
start to appear. 



V. STEEP EDGES 

There is a mystery in the wave packet evolution yet to 
be explained, that is, the appearance of steep edges. Our 
following analysis shows that the formation of steep edges 
can also be understood in terms of the BEC dynamics in 
the double-well system. For future convenience, we write 
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Figure 11: The evolution of relative phases for each pair of 
neighboring wells in the optical lattice. N = 2000. 
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Figure 12: The evolution of relative phases for each pair of 
neighboring wells in the optical lattice. N = 5000. 
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Figure 13: The evolution of relative populations s for each 
pair of neighboring wells in the optical lattice. N = 2000. 



which are derived from the Hamiltonian in Eq.(|5j). 

In Fie llll we have plotted how the relative phases 9 
evolve with time for N — 2000. In our calculation, the 
relative phase 9 is defined as the phase difference between 
the middle points in two neighboring wells. Initially, the 
relative phase is zero for every pair of double-wells in 
the optical lattice as indicated by a horizontal line. As 
the evolution goes on, the line of relative phase begins 
to incline with an increasing slope. The tendency stops 
when the two end points of the line reach ±7r/2, respec- 
tively. This is around t = 80ms, right when the steep 
edges appear. As shown in Fig^J the situation is simi- 
lar for N = 5000. We also observe from Figs Mil that 
the relative populations in the middle of the wave packet 
remain largely zero before the appearance of steep edges. 

According to Eq.JfjJ, at 9 ~ 7r/2 and s ~ 0, the tunnel- 
ing or transfer of atoms between these neighboring wells 
is the largest. This means that there are more atoms 
flowing into some particular wells than going out, thus 
generating steep edges. We notice that right after the 
appearance of the steep edges, both s and 9 become 
rather "random" . This is likely due to the complicated 
dynamics caused by steep edges. 
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Figure 14: The evolution ol relative populations s for each 
pair of neighboring wells in the optical lattice. N = 5000. 



VI. CONCLUSION 

With the GP equation, we have studied the wave 
packet dynamic of a BEC in a one-dimensional optical 
lattices. We find an intriguing self-trapping phenomenon 
in the expansion of the wave packet, agreeing with a 
recent experiment [(J. Moreover, we find that the self- 
trapping is only temporary and the wave packet contin- 
ues to grow at long evolution times that are beyond the 
current experiment 6] . The analysis of our numerical re- 



sults shows that the self-trapping in the optical lattice is 
closely related to the self-trapping found in the system of 
a BEC in a double-well potential. We also showed that 
the steep edges appearing the wave packet evolution do 
not necessarily lead to self-trapping and they can also be 
understood in terms of the dynamics in the double-well 
systems. 
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Appendix A: COMPUTATION OF £ IN OPTICAL 
LATTICES 



We choose one well and its neighbor as a double well 
trap. The GP equation is as Eq. |@J| except the potential 
is replaced by V(x) = V cos(x) for |x| < 27t and V(x) — 
V for |x| > 2ir. Then we write the wave function as a two- 
mode wave function: <f> = aui{x) + bu2(x) and let |a| 2 = 
N a /(N a +N b ) and \b\ 2 = N b /(N a +N b ) (N a is the particle 
number in well a and N b the particle number in well b). 
Plugging the double mode wave function into the GP 
equation with potential V(x) and using the tight-binding 
approximation, we obtain the effective Hamiltonian: 



H eff = -c/2(|a| 2 - |6| 2 ) 2 + v(a*b + ab*), 



(Al) 



where the parameter c = (N a + N b )Ng J \ui(x)\ 4 dx and 
v = N f(-±Vu{(x)Vu 2 (x) +V'(x)ul(x)u 2 (x))dx. The 
wave function u\(x) (or 112(2:)) is obtained as a ground 
state from Eq.Q with single well: V(x) — V cos(x) for 
< x < 2ir and V(x) = V for other x values. 
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